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Abstract: It is usually assumed that dark matter direct detection is sensitive to a large 
fraction of the dark matter (DM) velocity distribution. We propose an alternative form 
of dark matter-nucleus scattering which only probes a narrow range of DM velocities due 
to the existence of a resonance, a DM-nucleus bound state, in the scattering - resonant 
dark matter (rDM). The scattering cross section becomes highly element dependent, has 
increased modulation and as a result can explain the DAMA/LIBRA results whilst not being 
in conflict with other direct detection experiments. We describe a simple model that realizes 
the dynamics of rDM, where the DM is the neutral component of a fermionic weak triplet 
whose charged partners differ in mass by approximately 10 MeV. 



1. Introduction 



Although the fact that a large fraction of the matter in the universe is non-baryonic is beyond 
doubt, the exact form of the dark matter (DM) and the nature of the dark sector to which 
it belongs is still shrouded in mystery. The DM puzzle is under assault simultaneously on 
several fronts: direct detection experiments search for DM in our galactic halo through its 
recoil off nuclei in the apparatus, indirect detection experiments look for DM by searching 
for the standard model (SM) particles produced in DM-DM annihilation in the galactic halo, 
and collider experiments hope to produce DM in high energy collisions and infer its presence 
from large amounts of missing transverse energy in events. We concentrate here on direct 
detection experiments. 

The recent results from DAMA /LIBRA BJ have clear evidence of an annual modulation in 
the rate of nuclear recoils in their apparatus. Taken together with their previous DAMA/NAI 
results ]|], they observe this modulation at the 8.2 a confidence level (C.L.). Furthermore, 
the peak of this modulation is consistent with June 2 nd which is also the peak in the DM 
velocity relative to the Earth. The interpretation of the annual modulation as resulting 
from DM recoiling off nuclei in the DAMA detector is appealing. However, for the case of 
a conventional weakly interacting massive particle (WIMP) scattering off nuclei, explaining 
the DAMA modulation by DM recoils predicts many recoil events at other direct detection 
experiments, which are not observed. This apparent tension ||, ||, |j between DAMA and other 
experiments, such as CDMS §, KIMS §, XENON @, CRESST @, and ZEPLIN [0, 0, 
has lead to new DM explanations where the DM is not a conventional WIMP. 

Inelastic dark matter (iDM) proposes that the DM-nucleus interaction is inelastic 
in nature, the outgoing DM particle is actually an excited state heavier than the DM by 
0(100) keV. This changes the kinematics of DM scattering as the DM must have sufficient 
kinetic energy to be able to up-scatter, which favors the high velocity tail of the DM velocity 
distribution, and experiments that involve heavy elements. In particular, for an appropriately 
chosen splitting, there may be no events in CDMS but there will be events in DAMA and 
the rest, all of which involve elements at least as heavy as iodine. The lower cutoff on the 
DM velocity necessary for an inelastic scatter also means that the amount of modulation 
in the signal is enhanced, while at the same time the spectrum has a suppression at low 
recoil energy. These effects | Q , |14| , [l5| ] help to remove the tension between the DAMA result 
and the null results of the other experiments, for an analysis that reaches a slightly different 
conclusion see [16|. In addition, the small splitting in iDM means that the excited state may 
be very long lived, and the down-scatters that occur when it recoils against a nucleus lead 
to a novel signal, well outside the usual region of interest that may nonetheless be visible at 
direct detection experiments [17, [Df] . 

A recently proposed alternative to iDM is that the interaction between the DM and the 
nucleus is momentum dependent pC[ |. Such momentum dependence may come about 
from a form factor that is present in the coupling of the DM with the SM. The different DM 
experiments probe different ranges of momentum transfer and if the momentum dependence 
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is such that outside the range probed by DAMA the form factor is close to zero then DAMA 
would have increased sensitivity relative to the other experiments. However, there is still 
considerable overlap in momentum transfer between the various experiments and so it is not 
possible, for instance, to arrange for no events in CDMS and only events in DAMA [19]. Such 
a scenario prefers DM of mass around 50 GeV and requires either a non-Maxwellian halo [1^ 
or spin-dependent couplings p0[ . 

In both of these approaches the scattering cross section is taken to be velocity indepen- 
dent. For a given recoil energy the scattering takes place for all velocities of DM greater 
than some lower bound, v m i n , set by the recoil energy; in the case of iDM v m i n is higher 
than for elastic scattering due to the mass splitting. We investigate an alternative possibility, 
that the DM-nucleus scattering is velocity dependent, and in particular that the scattering 
takes place through the production of a resonance, the scattering cross section then has a 
Breit-Wigner form. The resonant form of the scattering, resonant dark matter (rDM), means 
that for given DM and target masses only a narrow range of velocities, around the resonance 
velocity, will actually scatter and leave a signal in the detectors. Since only a small portion 
of the whole DM velocity distribution is probed there is increased modulation of the recoil 
spectrum, despite the scattering being elastic in nature. As emphasized in [13] this helps to 
reconcile the DAMA results with the null experiments. In the concrete realization of rDM 
that we present here the resonance is a bound state of a nucleus and a charged partner of 
the DM, and the mass of the bound state is very sensitive to the particular target nuclei in 
the direct detection experiment. The velocity of DM in the Earth's frame is restricted to lie 
in the range km/s <> vdm ^ 800 km/s. If the corresponding resonance velocity does not 
lie below ~ 800 km/s there will be no resonant scattering, and the signal of DM recoils will 
be greatly suppressed. This high element sensitivity means that, unlike the two approaches 
outlined above, it is possible in rDM that iodine is the only element for which direct detection 
could occur. 

Before focusing on the particular form of the resonance as a DM-nucleus bound state, we 
first describe rDM in a model independent fashion, emphasizing the basic features should a 
resonance exist, in Section We go on to show how the two experiments involving iodine, 
DAMA and KIMS, can be explained simultaneously and discuss in general the allowed param- 
eter space. Then in Section ^, we introduce an explicit model, based on the work of Pospelov 
and Ritz pip , that realizes many of the necessary features. We demonstrate the sensitivity 
to the target element and discuss the various possible signatures and relevant constraints in 
Section ||. Finally, we conclude in Section |5[ 



2. Resonance Effects 



In this section, we consider other possible effects, which can dramatically change the results 
from traditional DM-nucleus elastic scattering cross section calculations. The detection rate 
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per unit detector mass at a DM direct detection experiment is given by |22[ 

d 3 vl^-a N F\E R ), (2.1) 



dR _ N T m N p x f ^^f{v,v E ) _ L ,>, 



dE R 2 n 2 Nx m y 

where tun ~ A mp is the nucleus mass with rap the proton mass and A the atomic number; 
F(Ep) is the nuclear form factor and accounts for the fact that the cross section drops as 
one moves away from zero momentum transfer; the two-parameter Fermi charge distribution 
is used to calculate F(Ep) throughout this paper [23]; Nt is the number of target nuclei per 



unit mass, given by Nt = Na/A with Avogadro's number, Na = 6.02 x 10 26 kg -1 ; is the 
cross section to scatter of a nucleus, and p^x 1S the reduced mass of the DM-nucleus system. 
The DM mass is m x and we take the local DM density to be p x = 0.3 GeV/cm 3 . The velocity 
of the dark matter onto the (Earth-borne) target is v. The Earth's velocity in the galactic 



frame, ve, is the sum of the Earth's motion around the Sun [22] and the Sun's motion in the 
galaxy p4|] . We assume the WIMP velocity distribution is Maxwell-Boltzmann with velocity 
dispersion v$ = 220 km/s. Thus, 

m **) = T-^m e-V**'?"* ■ (2.2) 

As a function of time in the galactic frame, the Earth's velocity is ve ~ 227+14.4 cos [2ir 
km/s, with T = 1 year and to is around June 2 nd . The DM velocity distribution is cut-off 
at the galactic escape velocity. Thus, the upper limit of the integration in (|2.1| ) is given by 
\v + Ve\ < v esc , and the lower limit, since we will consider elastic scatters, is given by 



lm N E R 



2 • (2-3) 



V 2 ^x 

The current allowed range for the galactic escape velocity |^5[ is 498 km/s < v esc < 608 
km/s. For concreteness we set v esc = 500 km/s. Increasing this value slightly increases our 
allowed parameter space, but the general features remain unchanged. Because of different 
energy detection efficiencies for different detectors, a quench factor /„ is introduced to relate 
the observed recoil energy, Ep, to the actual recoil energy Ep, Ep = Ep/ f q . This allows one 
to convert Eq. ( |2.1[ ) to the experimental differential spectrums as dR/dEp = 1/ f q dR/dEp. 
For example, we take the quench factor f q = 0.085 for the iodine element in the DAMA 
experiment. 

In the usual calculation the nuclear cross section is related to the nucleon scattering 
cross section, a p , by, 



(Zf p + (A- Z)f n f & 
f 2 u 2 



&N = 7^ ^~ a V ) (2.4) 



where are the coupling strengths of DM to protons and neutrons and \x nx is the DM- 
nucleon reduced mass. Here however, we wish to work explicitly with the nuclear scattering 
cross section, and leave relating it to the microscopic Lagrangian to later, section || In 
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Figure 1: DM velocity distribution after angular integration in the summer (red) and winter (blue) 
for the usual Maxwell-Boltzmann distribution (dashed) and for the case with a resonance (solid) at 
450 km/s with width 150 km/s, the escape velocity was taken to be 500 km/s. 

particular, in the usual approach a p is velocity, and element independent. We will see that 
in rDM these statements are no longer true. 

In rDM, the DM or its gauge partner, forms a short-lived bound state with the target 
nucleus. The mass of the bound state is denoted as m r . In this case the DM-nucleus elastic 
scattering cross section has a resonant structure. In the non-relativistic limit, one has s = 
(m x + itin) 2 + m x m^v 1 . For ^/s close to the resonance mass, a familiar formula is obtained, 



where E = ^fs is the center of mass energy; s x and sn are the spins of the dark matter 
and the target nucleus; J r is the total angular momentum of the resonant bound state. In 
the non-relativistic limit, the scattering process is dominated by the s-wave, so a selection 
rule, J r = s~^ + J/v, applies to the accessible bound state. T t ^, x n is the partial width of 
the boundstate decaying into \ P ms N an d is a function of the centre of mass energy. The 
total width Tt t may be larger than this width due to the existence of other decay modes, 
we will discuss this in more detail in Section ||. The momentum of the DM in the center of 
momentum frame is k = [iN\ v - Note that if there exists more than one resonance, the cross 
section is the sum over all resonances, each given by (p.5|). Since the DM is non-relativistic, 
we can rewrite the cross section as a resonance in velocity, 



on = 



2 J r + 1 7T ^ r ^ x N 



(2.5) 



{2s x + l){2s N + 1) k 2 (E - m r ) 2 + T 2 J4 ' 




(2.6) 
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Figure 2: The ratio of modulation as a function of the resonance velocity v r . We take this ratio as 
the summer event rate minus the winter rate relative to the summer rate. The iVron is chosen to be 
150 km/s and v esc = 500 km/s. The black (solid), blue (dashed) and red (dot-dashed) lines are for 
widths of 6 = 1, 10,50 km/s, respectively. 



Here, the normalization is 



2 J r + 1 4tt 2 5 2 rj_^y 

ao (2^ + 1X2^ + 1)^^ Tl t ' [ - (} 



and the resonance velocity and width are given by. 

2 2 (m r - m x - m N ) . 



2 (mr - m x - mjy) j4 = T^ 



For a narrow resonance the widths are constant, but for a wide resonance one must take into 
account their dependence on the velocity. In the narrow resonance case with 4Ci) r , Eq. (^ 
is well approximated by a delta function and is 

crjv = cr 5(v 2 - v 2 T ) . (2. 



For the existence of a resonance to have an observable effect the resonant velocity, v r , 
must be low enough that some DM particles have this velocity and the resonance must be 
narrower than the range of the DM velocity distribution, 5 < v esc . This second condition 
means that the resonance is relatively narrow, Ttot ^ 1 MeV. If these two conditions hold 
then the presence of the resonance effectively picks out only part of the Maxwell-Boltzmann 
distribution, Fig. [j]. In doing so the resonance increases the amount of modulation in the 
signal, since the signal now comes from a narrow range of velocities over which the summer 
and winter rates are discrepant, and there is no averaging over the whole velocity distribution. 

The total counting rate can be separated into two parts S = So + S m cos [2ir(t — to)/T], 
with Sq as the unmodulated rate and S m as the modulated rate. The amount of modulation, 
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Figure 3: Allowed regions of parameter space. The left plot has v r = 725 km/s and the right plot 
has m x = 500 GeV. In both cases the shaded contours denote the value of \ 2 f° r a fit to the first 12 
bins of the DAMA modulated data, and the region between the solid (dashed) black lines exceed the 
DAMA unmodulated (KIMS) data at 90% C.L. 



and the sign, now depend very sensitively on the values of v r and 5. This feature of increased 
modulation is similar to inelastic DM [12, 13], which also happens due to a restriction on 
which part of the Maxwell-Boltzmann distribution is probed, but occurs here despite the 
scattering being elastic in nature. 

In Fig. |^ we show a plot of the fraction of modulation, defined as the ratio of the difference 
of the summer and winter rates to the summer rate, as a function of the resonance velocity 
for three choices of the resonance width. For illustration purposes we take the DM mass to be 
500 GeV and we consider the ratio in the bin that corresponds to 3 keVee < Er < 3.5 keVee 
at DAMA. This is representative of what happens in other bins, and for other DM masses. 
For this bin, the corresponding f m i n is 150 km/s. 

From Fig. ^, we see that with sufficiently small 5 the ratio can be as large as 100%. If the 
resonance velocity is below the peak of the Maxwell-Boltzmann distribution it is possible to 
have a modulation ratio that is negative (meaning the winter time has more data than summer 
time) since below the peak the distribution is higher in the winter than in the summer. 

One candidate for the resonance is a DM-nucleus bound state, we will describe a model 
with just such an object in Section ||, in which case the resonant velocity may be highly 
element dependent. This opens up the possibility that iodine maybe the only element with 
an open resonance, and in other experiments the signal rate would be suppressed by ~ 
(5/v r ) 4 < 10 -4 . In this case only the DAMA and KIMS experiments, which both contain 
iodine, would be expected to observe a signal. We will show that because of the modulation 
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Figure 4: Left panel: Model predictions from a single resonance and the DAMA modulated date. 
Right panel: Same as the left panel, but for the unmodulated rate. For both plots, we have used 
m x = 500 GeV, <jq = 1-34 pb, v r — 725 km/s and 6 = 40 km/s. This model point is also allowed by 
the KIMS experiment. 



enhancement in rDM it is possible to simultaneously explain the modulated DAMA results 
without contradicting either the DAMA unmodulated results or the KIMS results. 

Independent of the specific details of a model, we have four parameters to describe the 
modulation spectrum . m x , v r , 5 and do- Because of the uncertainty in the number of 
parameters in a realistic model, we present the total x 2 to describe the goodness of fit for 
the first 12 bins, 2keVee < Er < 8 keVee, of the DAMA modulated spectrum. We also 
consider constraints from the "single-hit" unmodulated spectrum in DAMA and the results 
from KIMS experiment. The constraints from other experiments are model-dependent and 
will be discussed in Section |j| 

The allowed parameter space is shown in Fig. ||], where we treat the overall cross section 
o"o as a floating parameter to minimize the x 2 for a fit to the DAMA modulated spectrum, for 
a given set of m x , v r and 5. The 90% C.L. exclusion region, enclosed by the black solid and 
the black dashed lines, are from the unmodulated data of DAMA and KIMS, respectively. 
The first 22 bins (up to 8 MeV) of DAMA and all 8 bins of KIMS are included in this analysis. 
A constant quench factor, 0.085, for Iodine is used for the DAMA experiment. For KIMS, 
we interpolate the energy dependent quench factor [§§ as f q (E R ) = 0.1 e -° m35El{ + 0.06 for 



E R in keV, or equivalently, f q {E R ) = e -(E 2 R +5E R )/90 + Q Q6 for Er in ^eVee. We take the 
statistically averaged values of the four crystals in KIMS as the experimental input. 

From the left panel of Fig. ||, we can see that the fit to DAMA modulated data is almost 
independent of the DM mass, once m x is above 300 GeV. The KIMS constraint is stronger 
than the DAMA unmodulated one. Together, the parameter space with a total x 2 below 10 
is almost ruled out, but a large area of parameter space with x 2 between 10 and 15 is allowed. 
One can also see that a larger v r and a smaller <5 can more easily evade the constraints, which 
is due to the enhanced modulation effects for this part of parameter space. For illustration 



1 In the concrete model described in Sec. some of these parameters are related and there are only two 
independent parameters in total. 
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Figure 5: The same as Fig. |, but for m x = 500 GeV, a = 0.025 pb, v r = 470 km/s and 5 = 
0.075 km/s. 

purpose, we show one allowed point of the parameter space in Fig. ^, which has x 2 = 9-5, and 
another in Fig. ||, which has x 2 = 15.2 . For Fig. ||-[| we have used v esc = 500 km/s and, for 
DAMA, have taken into account the detector energy resolution, by smearing with a Gaussian 
distribution with a{E)/E = 0.U8/VE + 0.0091 with E in keVee. 

3. A Model 

As a realization of rDM we consider the case of a fermionic WIMP, x , which is nearly 
degenerate in mass with a charged partner, x > [ 21 1 , we will be interested in splittings of 



order 10 — 100 MeV. We also introduce a Z2 symmetry under which the dark matter particle is 
odd and thus x°, the lightest parity-odd particle, is a stable particle. In order to suppress the 
spin-independent coupling of %° to nuclei through Z boson exchange we take x to transform 
as (3,0) under SU(2) W x £7(l)y. 

The splitting of the charged from the neutral component can come from two sources. 
There may be higher dimension operators such as (xT a x)(H^T a H) suppressed by some scale 
A, which contribute v 2 /A = Ajjv- There are also loop generated contributions. After elec- 
troweak symmetry breaking the charged components of x are split from the neutral component 
by electromagnetic radiative corrections, the size of this splitting |28|, ^] is independent of 
the dark matter mass for m x 3> M\y and is around c^M^ sin 2 166 MeV. The higher 

dimension operator is comparable to the radiative correction for A ~ 10 6 — 10 7 GeV. For 
latter analysis, we require a mild cancellation between those two contributions and have 

m x ± - m x = A = 166 MeV - A uv ~ O(10) MeV . (3.1) 

The small splitting between x^ an d x° results in a long lifetime for the charged state. The 
width for the decay x^ ~^ x° e±iy is 

T m —^GlA 5 m ( ^ x 2 x 10~ 13 eV , (3.2) 
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neglecting the electron mass. The lifetime of x^ is t ~ 3 x 10~ 3 s for A = 15 MeV. Because 
of the almost degenerate masses of x anci X°> they have approximately the same amount of 
thermal relic abundance. The charged particles decay to x° shortly after the freeze-out time. 
To satisfy the observed value of dark matter relic abundance, the dark matter mass should 
be around 2 TeV [28]. Keeping in mind that the dark matter abundance may also originate 
from other non-thermal processes, we will treat the dark matter mass as a free parameter in 
this paper. We also note that the lifetime of x ± is far less than one second. Therefore, Big 
Bang nucleosynthesis is unaffected in the model. 

The lifetime of x ± is fairly long from the collider point of view. If produced in a collider, 
X s1 does not decay inside the detector and leaves a muon-like track but with a different mass. 
The existing searches for long-lived massive charged particles (CHAMPs) at LEP2 impose 



a bound on the mass, m x ± > 99.5 GeV |3J at 95% C.L.. Both DO |T| and CDF (3| 
at the Tevatron have also performed searches for CHAMPs. The current strongest bound 
is from CDF who searched for a single, isolated, weakly interacting CHAMP within the 
muon trigger acceptance. For the case at hand, the CDF bound constrains the sum of the 
production cross sections cr(x + X _ ) anci < J (x + X°)- The exclusion limit on the dark matter 



m x ± > 121 GeV at 95% C.L. 
For A > O(MeV) the kinetic energy of the DM is insufficient to allow an inelastic up 
scattering to x^- However, it is possible that for certain elements, the DM may form a short- 
lived electromagnetic bound state with the nucleus of an atom. The elastic scattering of DM 
off a nucleus is depicted in Fig. |6|. The intermediate boundstate may be in its ground- 
state or in an excited state, depending on which state is within the accessible range of dark 
matter kinetic energy. Interestingly, an 0(100) GeV WIMP has kinetic energy comparable 
to the typical energy splitting of the energy levels of nuclei. So, if the binding energy is large 
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Table 1: The electric binding energy for the ground state of the boundstate composed of x~ an( i 
z+iN, the nuclear mass difference Am = tha ^ — niA N and the neutron emission energy S n . The 
DM mass is taken to be 500 GeV. 



enough to compensate for the splitting A, one or two bound states may be probed in the 
dark matter scattering process. 

For x° scattering off element the bound state forms between x~ an d z+i^- ^or 
large Z the orbital radius is smaller than the Bohr radius of the nucleus and x~ is bound 



inside the nucleus. Using the two-parameter Fermi charge distribution 23] for the protons 
inside the nucleus we calculate the binding energy of this system, by solving the Schrodinger 
equation [34]. For the case where the DM mass is much larger than the nucleus mass, the 



electromagnetic binding energies between x~ an d z + iN, for the relevant elements, are shown 
in Table |]. The (x~ ,z+i -^0 binding energy is small in comparison to the total binding energy 
of the whole nucleus and we assume that the presence of x~ inside the nucleus does not have 
a large effect on the charge distribution or the spectrum of nuclear excited states. The excited 
states of the bound state are ~ 1 MeV above the ground state, whereas the excited states 



of the nucleus alone are ~ 500 keV for a fixed angular momentum [35, pq ]. In order to 
conserve angular momentum, the excited states of the nucleus alone should have angular 
momentum as Jjv(z+i,A)* = Jn(z,A)> Jn(z,A) ^ 1 f° r the fermionic dark matter considered 
in this paper. If the excited states of the bound state can be probed, the selection rule is 
I J n(z+i,A)* + L + 1/2] = | J n(z,A) + 1/2], with L as the orbital angular momentum between 
X~ and i +l N*. 

The resonance may form in its ground state, or an excited state, and the nucleus itself 
may also be in an excited state. Expressing the total excitation energy of the resonance above 
its ground state as u, the resonance mass is given by, 



m r = m v - + mA w + oj — Eh , (3.3) 

A Z+l 

Thus, the resonance speed squared is 



2{A + m A N -mA N + u-E b ) 
v 2 r = ^ - . (3.4) 

The mass difference between neighboring elements, due to nuclear binding energy, depends on 
the elements and isotopes in question [35, [36] and for the elements of interest lies in the range 
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MeV < rriA N — m,A N ^ 14 MeV. The nuclear mass differences for the corresponding atoms 

Z+1 JV r-, z 

are listed in Table [y. The mass differences for nuclei are offset from the mass differences for 
atoms by approximately 511 keV. The kinetic energy of the reduced mass system, 1/2 p x N v 2 
is around 200 keV. If the splitting, A, is comparable to the other scales in the numerator 
of (3.4), then for some element there will be a bound state which is accessible by the DM. 



Since the binding energy is dependent on Z, and therefore element, it is possible that some 
targets will have a resonantly enhanced elastic scatter and others will not, and instead have 
a rate suppressed by 5 /vf. We discuss this element dependence, how it may simultaneously 
explain results of DAMA and the null results of other direct detection experiments, and its 
potential signals in more detail in Section [|]. 

Having discussed the resonance speed, we turn to discuss the various ways it can decay 
and estimate its decay width 5. The resonance may simply fall apart back into the initial 
state, namely x° an d ^iV, this is the only possibility if the bound state is made in its ground 
state. The width for the resonance to decay to the DM plus the original nucleus can be 
extracted from the cross section for x~ z+i^ to scatter into x° z-N, 



*XN 



(TV 



f x - (r) Pp (r) 4 vr r 2 dr % = G| , (f) ^ ^ ^ 



2 



Here, p p (r) is the proton charge distribution with the normalization J p p (r) 4ir r 2 dr = Z; 
f x -(r), normalized as f f x ~( r ) 4vr r 2 dr = 1, is the wavefunction of X— with respect to the 
nucleus. For iodine, the falling-apart decay width is calculated to be 

x Hl : T r ^ xN ps 0.006 eV . (3.5) 

If the bound state forms in an excited state it has the possibility to de-excite by emission 
of photons and/or neutrons. The width for emission of photons depends on the multipole 
moment involved in the transition. For nuclear electric transitions of multipole moment L, 



the width is given by [37], 



87r(L + l) / 3 \\ 2L+1 ^ 2L 
L[(2L + 1)!!] 2 \L + 3 / 



and magnetic transitions are given by 

w-iB ("'-i^) 2 4 x (i^) 2q < +1 '"" 2 ' < 3 ' 7 ' 

where typically around 3, is the nuclear magnetic moment of the proton in nuclear magne- 
ton units. We have assumed that the radial nuclear wavefunction has a step function profile, 
going to zero beyond the nuclear radius, r^. For a given L, magnetic transitions are negligible 
relative to electric transitions. Taking the nuclear radius to be ~ 6 fm |^] and a typical 
transition energy of 100 keV, we find 

T 7 (£l) w 0.2 eV, T 7 (£2) ps 3.8 x 10~ 9 eV, 

r 7 (Ml) « 4.8 x 10" 4 eV, T 7 (M2) « 7.4 x 10" 11 eV. (3.8) 
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Similarly, it is possible that rather than the nucleus changing energy level, the transitions 
take place due to x~ changing energy level. We find from solving the Schrodinger for excited 
states of the bound state, with the charge distribution of the nucleus as before, that the 
typical radius of the x~ orbit is 1 fm, and the typical transition energy as u ~ 500 keV — 1 
MeV. We assume that there are no magnetic transitions, and the widths for photon emission 
due to electric transitions of x~ are comparable to those of the nucleus ( p.8| ). Note that if the 
emitted photon energy is smaller than the kinetic energy of the incoming DM, it is possible 
that the excited states emits a photon and subsequently falls apart into the x an d the original 
nucleus. 

For a heavy element the neutron and proton separation energies are comparable with 
S n ~ S p ~ 0(5 — 10) MeV. Due to the Coulomb barrier proton emission is greatly suppressed 
relative to neutron emission. If the DM-nucleus bound state forms with the nucleus in an 
excited state whose energy is above the neutron separation energy a neutron will be emitted 
from the bound state. The width for this process is large T n ~ 0.1 — 1 MeV [38|. The emission 
of the neutron may leave the nucleus in its ground state, but if not subsequent emission of 
photons will de-excite the nucleus. Unlike the case of photon emission, where depending on 
the transition the photon energy can be less than the initial DM kinetic energy, the energy 
carried away by the emitted neutron is large enough that the bound state is energetically 
forbidden from falling apart and releasing the DM. The difference in nuclear binding of %N 
and ^JV is 0(MeV) and so is sufficiently large that the DM remains bound to nucleus ^ 1 N 
after neutron emission. 

Which of the possibilities outlined above occurs depends on the details of the spectrum 
of nuclear excited states and DM-nucleus excited states, which are highly element dependent. 
This allows for the possibility that only certain elements are capable of observing rDM whereas 
for others DM scattering is highly suppressed. In particular, rDM allows the results of DAMA 
to be compatible with the null results from other direct detection experiments. 



4. Element Dependence and Potential Signals 

One of the interesting features of the model of rDM described above is that it naturally leads 
to a high sensitivity to the target nucleus involved in the DM scattering. Thus, it is possible 
that only for particular elements, and therefore particular experiments, there is an accessible 
bound state and that for others there is no rDM scattering. To make definite predictions 
requires detailed knowledge of the available states of the DM-nucleus bound state. However, 
even without this detailed knowledge there are several general statements that can be made. 

The experiments searching for direct detection of DM fall into two classes: those which 
contain only low-Z (Z< 40) elements - CDMS (Si, Ge) Q, and those that contain at least one 
high-Z element - DAMA (Na, I) Q, KIMS (Cs, I) §, XENON (Xe) [§, CRESST (W) §, 
and ZEPLIN (Xe) ]i~0| , O]. There is a large difference in binding energies between the low-Z 
and high-Z elements, see Table [j], so that if the splitting A is large enough, A > 11 MeV 2 , 

2 2.1 MeV is the mass difference of 74 As and l\Ge. 
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there will be no available bound state at CDMS but there may be bound states available 
at the other experiments. In this way the rate at CDMS is too small for any signal to be 
observed. 

Presumably iodine is an element with an available bound state, since DAMA has observed 
a modulated signal. If the resonance takes place in the ground state of the xenon-DM system, 
or if the excited state can not decay via an electric or magnetic dipole transition Eqs. ( |3.6| - 
p.7| ), so that the dominant width is for the bound state to fall apart into DM and the original 
nucleus. Then we estimate the resonance to be very narrow, we refer to this as Case I. 
Case I: 



tot ~ 1 r^x N 



0.006 eV . 



'"' - 0.075 km/s, 



2J r + 1 4vr 2 5 2 , , . 

= fn ,i /o ~T\ 2 ~ °- 025 P b > 41 

(2s x + 1)(2sjv + 1) IJ-xN v t 

for sn = 5/2, s x = 1/2, J r = 3, v r = 470 km/s and m x = 500 GeV. These are the numbers 
used to generate Fig. ||. In this case, the dark matter elastic scattering is the main process. 
The inelastic process with extra gamma rays emitting is suppressed by T^(E2) /T r _^ x ^ ~ 
10 -4 . We focus on this case as the explanation of the DAMA signal. There are only two model 
parameters m x and A (or v r ). So, we have x 2 /d.o.f = 15.2/10 = 1.52, which corresponds to 
a p- value p = 0.12. 

Alternatively, photon emission process may dominate, we refer to this as Case II. 
Case II: 

,„, I - % 0.2 x ( — — ) '' eV . 

7 V 100 keV 



0.1 x 1— ) km/s 



y fi xN V 100 keV 

2J r + l 4vr 2 ,5 2 r 2 N / a; x-3 

^0 = 77, ; TTTo ~~rT ~o a t^2 w L4 x 10 x mn i T7 P b - ( 4 - 2 ) 



{2s x + l)(2 SiV + 1) V 2 xN v 2 r T 2 V 100 keV 

for sn = 5/2, s x = 1/2, J r = 3. The DAMA modulated spectrum can also be fitted by 
choosing uj^ = 50 keV, v r = 450 km/s and m x = 250 GeV (x 2 = 15.6 and p = 0.11). 
If the photon emission width does not completely dominate the width to fall apart, this 
situation may still be capable of explaining the DAMA excess, and in addition there will be 
correlated signals which can be searched for. Since the elastic dark matter scattering is not 
the main process, Case II predicts photons with a few hundred keV energy at DAMA. These 
gamma rays should have energy matched to a nuclear spectral line of the element l2 \Xe. The 
number of gamma rays events will be larger than the number of modulated events by a factor 
of r 7 /r r ^ x 7\r ~ 10 — 100. If the energy of the emitted photon is always below the kinetic 
energy of the incoming DM the bound state can emit a photon and then fall apart. Again, 
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along with the DAM A signal we would expect photon lines. But if the emitted photon energy 
is larger than the DM kinetic energy then the DM will remain permanently bound to the 
nucleus, and it is not possible to explain the DAMA excess. 

The binding energy grows with Z and it is possible that for some high-Z element the 
difference — A is large enough that the resonance, if it occurs, is above the neutron 
separation energy S n ~ 5 — 10 MeV; the process of proton emission from the nucleus requires 
tunneling through the Coulomb barrier and is suppressed relative to neutron emission, we do 
not consider the case of proton emission here. For such an element a dark matter scatter has 
very different kinematics from the usual elastic scatter. The recoil energy is determined not 
only by the kinetic energy of the incoming dark matter ~ 100 keV but receives contributions 
from the energy released when the neutron is emitted, Ef, — A — S n ~ MeV. Such a recoil is 
outside of the usual energy range searched for at conventional direct detection experiments 
and may require dedicated experiments to find. If the resonance level is above the neutron 
separation energy this decay channel dominates and the DM remains bound to the nucleus. 

While it is possible, for 11 ^ A/MeV <; 16, that there will be no events at CDMS 3 
and the scattering off tungsten will involve neutron emission and likely be unobservable at 
CRESST. The proximity of binding energies for experiments involving iodine, caesium and 
xenon present more of a challenge. KIMS in particular contains the same element as DAMA 
and there is no confusion as to the position of the resonant velocity. As a result the KIMS 
results place strong constraints on rDM, see Fig. ||. 

The xenon based experiments of ZEPLIN and XENON may also present strong con- 
straints. The splitting of the appropriate energy levels in 55CS is ~ 500 keV, so it is possible 
that no resonance is available for these xenon experiments. However, it is difficult to be sure 
due to the comparable binding energies of iodine and xenon and the preponderance of different 
isotopes of xenon present in the detectors. For simplicity we have concentrated throughout 
on one of the high abundance isotopes. Furthermore, rDM has the feature of increased mod- 
ulation (see Fig. ^) and the XENON data was taken between October and February, which 
somewhat weakens their bounds on rDM. Similar to the case II for DAMA, the spectral lines 
corresponding to the isotopes of cesium may be the dominant signals at the ZEPLIN and 
XENON. A more detailed study of the spectra and number of the predicted events in other 
experiments and the allowed resonance speeds for each experiment is warranted. 

It is amusing to consider the possibility that the signal in DAMA is not coming from 
scattering off iodine but instead from scattering off sodium. For A < 2 MeV it is possible 
that in all other experiments, with the exception of silicon for which there is no resonance, 
that there is neutron emission. We have been unable to find parameters within our simple 
model that allow this to work whilst keeping the dark matter mass above the collider bounds 
mentioned earlier, but it is an intriguing possibility within the rDM framework. 

3 Although there is no tree-level elastic scattering there can be scattering at the loop level. The rate for this 
is highly suppressed, a v ~ 10~ 45 cm 2 , and will not be observable until the next generation of DM direct 
detection experiments. 
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4.1 Other Constraints and Signals 

The process of DM capture with neutron, or photon, emission may occur for many high-Z 
elements, and so searches for anomalously heavy elements place a constraint on models of 
this type. However, very few of the searches for anomalously heavy elements have been with 
high-Z elements, there have been some searches for heavy nuclei of gold and iron using mass 
spectroscopy [39] with bounds on the allowed fraction of fAu ^ 0(1CT 10 ) and fp e <^ £?(10 -8 ), 
for a 500 GeV DM particle. Assuming that an appropriate energy level exists the expected 
fraction of element X that will have captured a DM particle during exposure time r, taken 
to be of order the age of the Earth, is 

is/ (°v) \ ( t \ /500 GeV\ fi . 

where for the ambient dark matter density we take n x = 0.3 GeV/m x cm -3 . Using the 
delta-function limit of (|2.6|) we find that the effective cross section is given by, 



Over the whole parameter space shown in Fig. || the fraction of heavy elements that have 
captured DM is just below the present bound. For example, the Case I in Eq. ( |4.1| ) has 
/x~4x 10~ 12 . However, an improvement in the bound would start to probe much of the 
parameter space. Moreover, an extension of these searches to other high-Z elements would 
shed light on the splitting A and which elements have accessible energy levels. 

In the simple model described in Section || the DM is part of a weak triplet. As mentioned 
earlier the present collider bounds, from the Tevatron, on the charged state in the triplet puts 
the mass scale above ~ 120 GeV. For masses close to this bound, m x = 150 GeV, the cross 
section for pair production of the charged states at the LHC is large. At s 1 / 2 = 10 TeV it 
is 180 fb and at s 1 / 2 = 14 TeV it is 280 fb, but it drops to ~ 1 fb for m x = 500 GeV at 
s 1 / 2 = 14 TeV. It can be searched for as long-lived charged states at the LHC. 

The main dark matter annihilation products are gauge bosons, and the typical 
annihilation cross section is order of 1 pb. However, if the splitting, A, is small enough that 
there is an available (x + X~) resonance the DM annihilation cross section may be greatly 



enhanced 140], 41, J21fl. In this case, lots of positrons, anti-protons, photons and neutrinos can 



be produced and contribute to cosmic rays. 



5. Discussions and Conclusion 

Resonant DM provides a new scenario to explain the DAMA modulation results, while being 
consistent with other dark matter direct detection constraints. The model presented in this 
paper is very simple: a dark matter particle, which is part of an electroweak triplet which 
has an O(10) MeV mass splitting between its neutral and charged components, can realize 
all the features necessary for rDM. In our model the resonance present in the DM-nucleus 
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scattering is an electromagnetic bound state of the charged partner of the DM and a nucleus. 
The simple model presented here may not be unique, other models can also be constructed 
to realize rDM, if there contain a bound state of dark matter and the nucleus. 

There are two novel features of rDM. Firstly, the resonance effect picks out from the 
Maxwell-Boltzmann dark matter velocity distribution a narrow window of DM velocities 
around the specific resonance velocity. If this velocity happens to be in the high-speed tail 
of the velocity distribution, the modulation amplitude of DM scattering can be enhanced. 
Secondly, rDM depends on many detailed properties of nuclei. Other than the mass, Z and 
A of a nucleus, the resonance speed can also depend on the energies of its various excitation 
levels. It is conceivable that the only element capable of direct detection of the DM in the 
halo is iodine. Because of this the KIMS experiment is in an ideal situation to confirm or 
deny the DAMA excess and the rDM model. 

The elemental dependence of the scattering rate in rDM means that it is not straight- 
forward to make precise predictions for the event rate at other experiments, it would be 
worthwhile to address the allowed parameter space for the other experiments in a model in- 
dependent fashion. However, it is straightforward, to accommodate the null results at CDMS. 
The model predictions for ZEPLIN and XENON depend highly on a more precise calculation 
of the binding energy and more accurate values for nuclei energy levels. For CRESST, de- 
pending on the mass difference of the DM particle and its partner, one can either have lots of 
MeV neutron events or a suppressed number of nuclear recoil events if there is no rDM effect. 
Exactly what occurs at DAMA and KIMS depends on which state the rDM effects probe. 
We have two cases: case one has only elastic scattering, if the electric and magnetic dipole 
radiation is forbidden, for instance if the bound state forms in its ground state. While case 
two predicts lots of gamma rays in addition to the elastic scattering events. The energy of 
those gamma rays should correspond to a nuclear transition in l J[Xe. It would be interesting 
to search for this signal at existing direct detection experiments. 

We have surveyed the allowed parameter space of rDM in a model independent fashion 
and have found it compatible with the DAMA and KIMS data for a wide range of DM 
mass, m x > 100 GeV. We have focused on the case of one resonance but note that there 
can in principle be several resonances available. For our concrete model, where the bound 
state requires a charged state nearly degenerate with the DM, there are colliders bounds on 
the charged state requiring m x ^ 120 GeV. This charged partner can be searched for at 
the Tevatron and the LHC. In addition, the DM can form stable bound states with heavy 
elements, the capture rate is just beyond the present bound. However, this search has only 
been carried out for a few heavy elements and by extending these searches for anomalously 
heavy nuclei to a larger set of heavy elements one can place strong constraints on rDM 
and probe the element dependence. In rDM models, the search for heavy nuclei is highly 
correlated with DM direct detection searches. Discovery of an anomalously heavy element 
would identify what material should be used for future direct detection experiments. 

In conclusion, the resonance effects can dramatically change the "traditional" dark matter 
elastic direct detection calculation. Resonant DM can enhance the modulation effect and 
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relies strongly on the detailed nuclear properties of different elements. The modulated data 
at DAMA can be explained and the constraints from unmodulated data can be satisfied. 
Furthermore, the rDM can reconcile the apparent contradiction between DAMA and other 
experiments like CDMS, XENON, KIMS, ZEPLIN and CRESST. We have described a simple 
model, which can realize all the features of rDM. It has dark matter, with m x > 120 GeV, to 
be part of an electroweak triplet which has a mass splitting between its neutral and charged 
parts of around 10 MeV. 
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